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Abstract. A model is presented for the chemical evolution of the solar neighbourhood 
which takes into account three families of galactic objects, according to their condensation 
states: stars, refuses and gas. Stars are defined as every condensed objects with masses 
greater than or equal to the minimum mass which ignites hydrogen and which will give rise 
to an evolutionary track on the HR diagram to the left of Hayashi's limit; refuses include 
the remnants, which are compact objects resulting from stellar deaths, and the residues, 
which have masses not large enough to ignite hydrogen; gas is defined as the mass which 
can be condensed to form stars and/or residues. We have developed equations for the 
mass evolution of each family, and have studied the gas metallicity distribution within 
the framework of the instantaneous recycling approximation, adopting different initial 
conditions. In order to constrain the model parameters we have also used preliminary 
evaluations of comet cloud masses to investigate the role of the residues as sinks of heavy 
elements in the Galaxy. 
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1. Introduction 



Models for the chemical evolution of galaxies usually include only two classes of objects, 
namely stars and gas (cf. Tinsley, 1980). This classification has as main argument the 
simplification of the equations, and generally does not cause significant difficulties in the 
derivation and interpretation of the quantities which are effectively compared with the 
observational data. 

On the other hand, more complete formulations already appeared in some of the 
early works on chemical evolution (Schmidt, 1959) and also in more recent treatments of 
the evolution of our Galaxy (Tinsley, 1981, Rana and Wilkinson, 1986) and other spiral 
galaxies (cf. Ferrini et al., 1992). 

The main characteristic of these treatments is the inclusion of non-stellar objects, as 
our present knowledge makes us sure that an important quantity of metals is probably 
locked up in some galactic objects such as planets, comets, etc. (cf. Bailey, 1988). Comets 
are particularly interesting in this respect, as they seem to be a very common phenomenon 
associated with star formation out of a protostellar cloud (Vanysek, 1987a,b). Similar to 
the case of interstellar grains, comet formation has probably a small effect in the mass 
balance of the Galaxy (cf. Meusinger, 1992). On the other hand, these objects may affect 
some observational properties such as the extinction and polarization of visible light (cf. 
Greenberg, 1974), and can in principle have an influence on the chemical evolution of the 
solar neighbourhood as metal sinks, as suggested by Tinsley and Cameron (1974), Vanysek 
(1987a,b), and Stern and Shull (1990). 

In the present paper, we have introduced a consistent treatment of the evolution- 
ary histories of the different families of galactic objects, taking into account the following 
condensation states: stars, refuses, and gas. We then study the derived metallicity dis- 
tribution for the one-zone model of the solar neighbourhood considering a set of initial 
conditions, and compare our results with observed data from stars. Finally, we make some 
preliminary calculations of comet cloud masses, in order to investigate the role of comets 
as heavy element sinks in the galactic disk. 

2. Basic Equations 

The adopted families of stars, refuses and gas are associated with the following classes of 
objects, respectively: 
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Stars, which are defined as every condensed object formed with masses m > mj, where m; 
is the lower limit for the stellar masses, or the minimum mass which produces hydrogen 
ignition and which will give rise to a track on the HR diagram to the left side of Hayashi's 
limit. 

Refuses, which include remnants, or compact objects resulting from stellar deaths, and 
residues of star formation, which are objects condensed from the gas, with masses in the 
interval m g < m < mi, where m g is the maximum assumed mass for the gas (e.g. grains) . 

Gas, which is the mass that can be condensed to form stars and/or residues. 

Adopting a model for the chemical evolution of the solar neighbourhood with no infall, 
the total mass of the system is constant and given by 

M = M g + M s + M r (1) 

where M g ,M s , and M r are the total masses in gaseous, stellar and refuse condensation 
states, respectively. The gas (/z) and refuse (k) mass fractions are defined by 



and 



so that 

M 8 = (l-n- k)M. (4) 

We will adopt the usual sudden mass loss approximation, where the stars undergo the 
entire process of mass loss after a well-defined lifetime. Let w m be the remnant mass of a 
star with initial mass m and lifetime r m . The rate of mass locked up in the remnants, due 
to the death of the stars which were born at instants given by (t — r m ) is obtained by 

L(t)= w m {m)y(t-T m )<$>(m)dm (5) 

J m t 

and the total ejection rate due the death of these stars is 

E{t)= [m-w m {m)]^{t-T m )^{m)dm, (6) 

J m t 
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where m t is an appropriately chosen turnoff mass; m u is the upper limit to the stellar 
mass, or the maximum mass admitted to stars; \1/ is a generalized formation rate, defined 
as the total mass condensed into galactic objects per unit time, and $ is a generalized 
initial mass function, normalized as 



where x\ = 0, = m g , X3 = m/, £4 = m t , x$ = m u , and p, 7, ( and e are, respectively, 
the first, the second, the third and the fourth terms in the sum. Each of these four terms, 
multiplied by the generalized formation rate will give respectively: the formation rate of 
objects with masses m < m g , considered as gas, p^(t); the formation rate of the residues, 
7*(t); and the stellar formation rates, C^(t) an d s^(t), relative to stars with masses in 
the intervals mi < m < mt and mt < m < m u , respectively. 

It is worth noting that the generalized initial mass function is not necessarily continu- 
ous in the above intervals. However, we have assumed its continuity and will take m g ~ 0, 
so that p ~ 0. We will assume that every residue is gravitationally tied with a star and 
that the formation of stars with mass m > mt is not accompanied by formation of residues 
(Stern and Shull, 1990). 

Some residues may undergo evaporation of their H and He. If we assume that 7 is 
the fraction of the generalized formation rate that will initially produce such residues, the 
gas will be replenished by a mass per unit time (1 — Z)^ due the evaporation of H and 
He, where Z is the metallicity of the gas and Z^ is the mass of metals which go into 
these residues per unit time. We will assume an instantaneous evaporation. Based on 
Tinsley and Cameron (1974) and Vanysek (1987), we set 7 > 7 - 7, so that as a first 
approximation 7^7. 

Adopting the instantaneous recycling approximation (IRA), r m 0. The ejection rate 
to the interstellar medium and the rate of mass locked up in remnants can be simplified as 





E(t) = m(t) 



(8) 



where R is the returned fraction to the interstellar gas, 



R = / [m — w m (m)]$>(m)d, 



m 



(9) 



and 



L(t) = (e-R)V(i). 



(10) 
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Recalling that p m in equation (7), that 7 ~ 7, and assuming further that Z C 1, we 
can write for the mass rates 

l Mg (t) = -(C + e - RMt) (11a) 

j t M r {t) = (e - R)*(t) (11c) 

Equation (lib) can also be obtained assuming a constant formation rate. In the framework 
of the IRA, it can be easily interpreted: the rate of change of total stellar mass is only due 
to the stars that live forever, which have masses in the range mi < m < m t . 



3. The generalized formation rate and initial mass function 

The slow rate of growth of the abundances of the heavy elements produced by the metal 
sink effect due to the refuses (Tinsley and Cameron 1974; Vanysek 1987a,b) is built in our 
model, and can be obtained by appropriately chosen fractions of the generalized formation 
rate. In order to determine this rate, we will follow Tinsley and Cameron (1974) and 
assume that the mass of metals which go into comets per unit time is at least equal in 
magnitude to the mass of metals which go into the associated star. Since we have assumed 
7 ~ 7, it follows that 7<, In order to estimate the fractions of the formation rate, we 
have used the stellar IMF from Miller and Scalo (1979) for m > m\. We have assumed 
that the generalized initial mass function ($) for residues is proportional to m~ x , and 
re-normalized $ in the interval (m g ,m u ) assuming 7< (. We have taken mi = O.1M 
(Larson, 1992), m t = 1M© and m u = 100M© (Tinsley, 1980). It can be shown that if 
7< C, then the slope x of the generalized initial mass function for residues should be < 1.8. 

The fractions of the generalized formation rate can then be computed and we have 
obtained 7 « 0.3, and Q « 0.3, so that e « 0.4. The returned fraction can be computed 
from equation (9) as R 0.24, where we have taken w m = 0.7 M Q for m < 4M , and 
w m = 1.4M© for m > 4M (Tinsley, 1980). 
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4. Metallicity Distribution 



Following Tinsley (1980), the metallicity of the gas for the conditions adopted here can be 
obtained from the equation 

j t Z(t)M g (t) = -Z(l -R)y + y ( £ + (-R)y (12) 

where we have kept the definition of the heavy element yield y as the total mass of new 
ejected metals relative to the mass locked up in stars and remnants, 

y=— ^ mpzm(m)$(m)dm (13) 

C + e - R J mt 

where pz-m is the so-called stellar yield, or the mass fraction of a star with mass m that is 
converted to metals and ejected. 

From equations (11a) and (12) the metallicity can be integrated as 
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+ Z ( A ) . (14) 



where \i = fj,(t), so that \x = hq and Z = Z for t = 0. We set t = in that instant when 
the disk reached its final mass M. 

In order to show the dependence of the heavy element abundance with 7, we will take 
y = 0.01 (Tinsley and Cameron, 1974; Maciel, 1992) and R = 0.24 as discussed in section 
3. We have also adopted an initially unenriched gas, so that Z w 0. Figure 1 shows Z as 
a function of the ratio fi/fio for some representative values of 7. We see that the smaller 
is the value for 7, the greater is Z{ji — > 0). Tentative limits for 7 are provided, assuming 
that Z{ji — > 0) w Z Q = 0.02 ± 0.01 in equation (14). On the basis of the assumed error, 
we can see in figure 2 that values for 7 in the range 0.19-0.38 are preferred. It can be seen 
that the residue mass fraction obtained in section 2 lies approximately in the middle of 
this interval. 

An analytical expression for the cumulative distribution of stars of a given metallicity 
can be derived for the one-zone model with metal retention by refuses. Recalling the 
definition of the gas and refuse fractionary superficial densities, equations (2) and (3), 
respectively, the fraction of stars born until the metallicity has reached a value Z is 

s (Z ) = = izm^m, (15) 
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where the subscript 1 indicates present values. From equations (11a) and (11c) we can 
write 

K = K + - -(jio-fl), (16) 

i — 7 — it 

where the subscript again indicates initial values. From equations (14), (15) and (16), 
and recalling the definitions (2) and (3) , we obtain, after some algebra 



a — b/j,o 
S(Z) = 



jz/Zi - 1) - {z/z 1 -Zp /Zj ) (jm /no ) t^t- 
z /z 1 -\ 



a — b\x\ 

where a and b are constants given by 



(17) 



a = 1 - k - - -no (18) 

l — 7 — K 

b -r4^R- (19) 

Of course, S is normalized, so that form (17) we have S = 1 for Z = Z\. 

In order to analyze the results of equation (17), we have varied the initial conditions, 
namely no, Ko e Zq. We have chosen eight sets of initial conditions, labeled by letters 
A to H, which are shown in table 1. To obtain numerical estimates, we have used the 
mass fractions given in section 3, namely 7 = 0.3, Q = 0.3, and e = 0.4. We have further 
assumed m = 0.1 (Tinsley, 1980; Pagel and Patchett, 1975), and R = 0.24. 

Figure 3 shows cases A to D. As a comparison with the observational data, the asterisks 
in the figure are obtained from the differential metallicity distribution of 132 G dwarfs in 
a cylinder passing through the Sun and perpendicular to the galactic plane (Pagel, 1989). 
We have taken Z\ = 1.19-Z©, according to Table 2 of Pagel (1989), which corresponds to 
the central value of the largest metallicity bin. 

Case A will arise in a disk with non unitary gas fraction, but with no initial refuses. 
The distribution will predict even greater values than the simple model (Schmidt, 1963). 
This is due to the fact that when we set no 7^ 1 and «o = 0, we are necessarily accepting 
some primordial stars with low metallicities. Case B is the simple model with metal sink 
effect. Cases C and D shows "prompt initial enrichment" models in which the burst of 
star formation will also lead to the formation of refuses. 

Figure 4 shows the effect of increasing Zq fixing the initial gas and refuse mass frac- 
tions. We can see that the fits to the observational data are much better, especially for the 
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models with higher initial heavy element abundances (models G, H). This result is partic- 
ulary interesting when we compare models A and E, where the inclusion of refuses and a 
very small initial heavy element abundance produce a large difference in the cumulative 
distribution at low metallicities. 



5. Comets and residues 

Tinsley (1974) has pointed out that two empirical results provide powerful constraints on 
chemical evolution models, namely the G-dwarf problem and the slow enrichment rate of 
the ISM. In our model, the assumption that comets are like sinks of metals explains easily 
this slow enrichment, provided that the slope of the generalized initial mass function for the 
residues is < 1.8. On the other hand, the G-dwarf problem is also explained if we postulate 
that a generation of primordial massive stars will give to the disk an initial metallicity, as 
well as some initial remnants. Of course, the initial conditions are connected, so that the 
values for (xq and kq are likely to depend on the value for Zq. The crucial assumption in 
our model is that 7^7. 

In order to obtain an order-of-magnitude estimate of the importance of comets as part 
of the residues, we have estimated the total initial residue mass 7W. Using 0.19 < 7 < 0.38, 
and a present value of the generalized formation rate similar to the star formation rate, \& ~ 
1OM pc" 2 Gyr" 1 (Tinsley, 1980; Miller and Scalo, 1979), we have 7^ ~ 1.9-3.8M pc" 2 
Gyr -1 . The corresponding term for comets can be estimated by 

Z ^ ~ N - h - M - (20) 
Vt 

where N c is the number of comets, h c is the comet galactic scale height, M c is the average 
nuclear mass of a comet, V is the total volume considered and r is the system lifetime. We 
have first taken into account the whole solar system, where N c ~ 2.5 x 10 6 (Allen, 1973), 
V ~ 3.1 x 10 -11 pc 3 and h c ~ 4.0 x 10~ 4 pc with a radius of 40 AU (Allen, 1973), r ~ 5 
Gyr, and M c ~ 1O" 16 M (Vanysek, 1987), so that ~ 6.5 x 1O" 4 M pc" 2 Gyr" 1 . 

Taking the average heavy element abundance during the solar system lifetime Z ~ 0.01, 
we have 7^ ~ O.O6M pc -2 Gyr -1 , which is much lower than the 7* fraction estimated 
above. Assuming now the existence of the so-called "Massive Oort Cloud" with h c ~ 0.10 
pc for an adopted radius of 10 4 AU (Vanysek 1987), V ~ 4.8 x 10" 4 pc 3 , r ~ 5 Gyr, 
and M c ~ 1O~ 16 M , we need N c ~ 10 13 comets to account for most of the residue mass, 
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in agreement with independent results by Stern and Shull (1990), Vanysek (1987a,b) and 
Greenberg (1974). 
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